Chaos-assisted directional light emission from microcavity lasers 
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We study the effect of dynamical tunneling on emission from ray-chaotic microcavities by in- 
troducing a suitably designed deformed disk cavity. We focus on its high quality factor modes 
strongly localized along a stable periodic ray orbit confined by total internal reflection. It is shown 
that dominant emission originates from the tunneling from the periodic ray orbit to chaotic ones; 
the latter eventually escape from the cavity refractively, resulting in directional emission that is 
unexpected from the geometry of the periodic orbit, but fully explained by unstable manifolds of 
chaotic ray dynamics. Experimentally performing selective excitation of those modes, we succeeded 
in observing the directional emission in good agreement with theoretical prediction. This provides 
decisive experimental evidence of dynamical tunneling in a ray-chaotic microcavity. 
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Dynamical billiards have served as a simple and 
generic model in studies of classical and quantum chaos. 
Whereas originally a dynamical billiard was introduced 
as an abstract model for studying the ergodic hypothesis 
nowadays it has been realized as a physical system 
that confines microwaves, acoustics, electrons and light 
0. As for billiards for light |^|, or optical microcavi- 
ties, the openness of a system intrinsic to light confine- 
ment necessitated the studies of decaying modes, or res- 
onances of the Helmholtz equation, which had been less 
studied in the field of quantum chaos. Resonance char- 
acteristics such as quality factors and emission patterns 
have been intensively studied in terms of classical and 
quantum chaos theory [6[. These studies are also moti- 
vated by practical applications of optical microcavities: 
When they are applied for laser resonators, there is a 
demand that modes have both high quality factors and 
high emission directionality. 

A high quality factor mode can be associated with a ray 
orbit confined by total internal reflection (TIR), whose 
representative example is the whispering-gallery (WG) 
modes of a disk cavity. It was demonstrated that the in- 
troduction of a slight asymmetric deformation to a disk 
shape enables the existence of a WG-like mode with emis- 
sion directionality [3]. However, in general, such a mode 
is embedded in a variety of high quality factor modes that 
do not necessarily have similar emission directionality, 
implying that experimentally it is not easy to selectively 
excite only desired directional modes. 

For a large deformation, ray dynamics inside the cav- 
ity becomes mostly chaotic, which causes degradation of 
quality factors. However, low- loss modes start to exhibit 
a universal emission pattern, which can be well explained 



by the unstable manifolds of an unstable periodic orbit 
close to the TIR condition [7]. The universal emission 
pattern can be highly directional by properly design- 
ing an unstable manifold structure [8[. Because low- loss 
modes, likely to be excited with a lasing medium, have a 
similar emission pattern, one can experimentally achieve 
directional emission without any settings for selective ex- 
citation. For those low-loss modes, relatively high quality 
factors are attributed to wave localization along TIR- 
confined unstable periodic orbits [8[, which is, however, 
generally weaker than localization along a TIR-confined 
WG orbit or a TIR-confined stable periodic orbit. 

In this Letter, introducing a novel largely deformed 
disk cavity, we study high quality factor modes with emis- 
sion directionality made possible by a dynamical tunnel- 
ing phenomenon [9|, LLOj . The key point is mostly chaotic 
ray dynamics having a single dominant stable periodic or- 
bit confined by TIR. The high quality factors come from 
strong localization of the modes along the TIR-confined 
stable periodic orbit. If one adopts a one-dimensional or 
integrable picture based, for instance, on the Gaussian- 
optic theory , one would expect extremely weak emis- 
sion due to TIR for those modes. However, in reality, 
tunneling from the periodic orbit to neighboring chaotic 
orbits takes place, where the latter eventually escape 
out from the cavity refractively following unstable man- 
ifolds. Namely, light is emitted with the assistance of 
chaos. Fabricating laser diodes with this deformed disk 
shape with an electrode contact patterned along the TIR- 
confined periodic orbit, we experimentally perform se- 
lective excitation of the high quality factor modes and 
demonstrate chaos-assisted directional emission in good 
agreement with theoretical prediction. 
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We first define the deformed disk cavity. In the polar 
coordinates, its boundary is defined by r(<p) = R(l + 
a cos 2(j) + b cos 4</> + c cos 60) , where is the size param- 
eter and a, b and c are deformation parameters fixed as 
a = 0.1, 6 = 0.01 and c = 0.012 [Fig. ffl[a)]. Features of 
internal ray dynamics can be best described by using the 
Poincare surface of section (SOS). The Poincare SOS is 
a reduction of ray dynamics to a two-dimensional map- 
ping describing successive bounces of a ray orbit. The 
variables of the mapping are a canonical conjugate pair 
called the Birkhoff coordinates, s and sin#, where s is 
the arclength measured along the cavity boundary and 
sin# the tangential momentum of an incident ray [Fig. 
Ufa)]. We show the SOS of our cavity in Fig. HJc), which 
reveals that ray dynamics is mostly chaotic, except for 
several islands of stability. Dominant islands are period- 
two islands at around sin = and period- four islands at 
around sin# = 0.7. FigureQJb) shows a rectangular orbit 
corresponding to the period-four stable periodic points. 
This orbit is confined by TIR when the laser is made from 
GaAs/AlGaAs quantum- well structure (i.e., the effective 
refractive index n = 3.3). In Fig. ffl^c), the critical line 
for TIR, sin# = 1/n, is plotted by a dashed line, above 
which rays are reflected by TIR. The key feature of our 
cavity is that there exists only one dominant stable pe- 
riodic orbit above the critical line and it is surrounded 
by a chaotic sea extending into the non-TIR (i.e., leaky) 
region, |sin#| < 1/n. 

We are interested in resonant modes corresponding to 
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the TIR-confined rectangular orbit. For those modes, 
there are two expected emission mechanisms. One is 
evanescent emission, which occurs at the four bounc- 
ing points of the rectangular orbit, yielding emission 
in the tangential direction to the cavity boundary, i.e., 
(j) = 45°, 135°, 225°, and 315°. Another mechanism is due 
to dynamical tunneling @-ll|. An orbit on the period- 
four islands tunnels to a neighboring chaotic orbit and 
it eventually diffuses into the leaky region, yielding re- 
fractive light emission. In order to capture the effect of 
dynamical tunneling on emission patterns precisely, next 
we study wave functions of resonant modes. 

Resonant modes can be obtained by solving a linear 
Helmholtz equation numerically [12| . We consider trans- 
verse electric (TE) polarization, in order to be consistent 
with our experiments described below. We systemati- 
cally study resonant modes whose scaled wave number 
kR is around 50, where k is the free-space wave number. 
We checked that properties of resonant modes discussed 
here are essentially the same for a larger cavity (e.g., 
fci?=100). 

Employing the Husimi projection of a wave function 
onto the SOS, we could classify resonant modes into 
three classes: rectangle-orbit modes, WG-type modes 
and chaotic modes. Here we exclude highly lossy modes 
such as the one residing on the period-two islands in the 
leaky region. Rectangle-orbit modes are strongly local- 
ized on the period-four islands, while WG-type modes are 
strongly localized on the upper or lower verges of the SOS 
(i.e., |sin0| £ 1). Chaotic modes are mostly distributed 
in the chaotic sea, even largely spreading into the leaky 
region. The first two classes of modes have high quality 
factors Q « 10 5 - 10 6 , as they are localized in the TIR 
regions, whereas typical chaotic modes have low quality 
factors Q « 100. 

We plot in Fig. [2]the wave function of a rectangle-orbit 
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FIG. 1: (color online) (a) Geometry of the deformed disk 
cavity, (b) The rectangular ray orbit, (c) The upper half 
of the Poincare surface of section for internal ray dynamics. 
The variable s is normalized by the total circumference. The 
four filled rectangles (■) plotted at sin# w 0.7 correspond to 
the rectangular ray orbit. The critical line for total internal 
reflection, sin# = 1/3.3, is plotted by a dashed line. 




FIG. 2: (color online) Wave function for a rectangle-orbit 
mode with a wave number KekR — 49.94 and Im kR — 
—0.00012. The intensity increases as the color changes from 
black to white. Outside the cavity, the intensity is plotted 
in log scale. The outline of the electrode contact window is 
indicated by white curves. 
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FIG. 3: (color online) The Husimi projection of the rectangle- 
orbit mode shown in Fig. [2] The intensity increases as the 
color changes from black to white. The intensity is plot- 
ted in log scale. The critical line for total internal reflec- 
tion is plotted by a solid line. Unstable manifolds emanating 
from a period-three unstable periodic point at (s,sin0) = 
(0.5,0.389...), indicated by a filled circle, are plotted by a 
solid curve. The high intensity spots around the critical line 
at s = 0.04 and s = 0.54 correspond to emission towards 
cj) — 90° and cj) — —90°, respectively. 



mode, where the intensity outside the cavity is plotted 
in log scale in order to clearly see the emission pattern. 
Figure [2] reveals that light is emitted towards (j) = ±90° 
on both sides of the cavity, i.e.,0 ~ 0° and 180°. By 
careful inspection, one finds that light is emitted from 
four cavity boundary points corresponding to polar an- 
gles (j) = ±4° and (j) = 180 ±4°. In the corresponding 
far-field pattern (to be presented in Fig. [5j), two distinct 
peaks appear at (j) = ±90° with a divergence angle 30°. 
Moreover, in each peak, we observe the interference pat- 
tern of light beams emitted from both sides of the cavity, 
whose oscillation period is given by A6 = 360° jkd with 
d = 2R (1 + a + b + c) being the width of the cavity. We 
confirmed that emission patterns are almost similar for 
all rectangle-orbit modes. Although it appears counter- 
intuitive that light is emitted from positions apart from 
the rectangular orbit, its mechanism can be explained by 
dynamical tunneling as we see next. 

Figure [3] shows the log plot of the Husimi projection of 
the wave function shown in Fig. [2j Because the Husimi 
projection is symmetric with respect to the transforma- 
tion (s, sin#) (1 — s, — sin#), only the upper half part 
is shown. The intensity is mostly concentrated around 
the period- four islands, but we can see some intensity 
also spread in the chaotic sea, which can be interpreted 
as the effect of dynamical tunneling. What determines 
the emission pattern is the intensity distribution in the 
leaky region. Around the critical line for TIR, we can ob- 
serve two relatively high intensity spots at s = 0.04 and 
s = 0.54. It is these that are responsible for the emission 
towards <j> ±90°. 

It is explainable by ray dynamics why emission occurs 
at s = 0.04 and s = 0.54. In Fig. [3j we superimpose 
unstable manifolds of a period-three unstable periodic 
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FIG. 4: (a) Scanning electron microscope image of a fabri- 
cated laser diode with the margin parameter M = 0.95 (see 
text), (b) CCD photo of a lasing diode with M = 0.95. (c) 
The wave function of a rectangle-orbit mode covered by the 
electrode metal area with M = 0.8. (d) CCD photo of a las- 
ing diode with M = 0.8. In (b) and (d), the photo is taken 
from a slightly oblique angle and the shade of an electrode is 
indicated by dotted lines. 



point onto the Husimi projection. This periodic point 
is located close to the critical line and thus its unstable 
manifolds govern a ray-dynamical flow around the critical 
line . The structure of the unstable manifolds indicate 
that ray orbits are transported from the TIR region to 
the leaky region at around s = 0.5 (and s = as well 
because of the 180° rotational symmetry of the cavity). 
This ray-dynamical argument also explains why emission 
patterns of the rectangle-orbit modes are almost similar. 

Now we would like to examine whether the chaos- 
assisted light emission can be experimentally observed. 
We fabricated single-quantum-well laser diodes with 
the deformed disk cavities discussed above. The laser 
diodes are fabricated by applying a reactive-ion-etching 
technique to a graded-index separate-confinement- 
heterostructure single-quantum- well GaAs/AlGaAs 
structure that was grown by metal organic chemical 
vapor deposition. The lasing wavelength is around 850 
nm. The cavity size R is set as 50 /im, yielding the 
scaled wave number kR = 370. In order to inject current 
only to rectangle-orbit modes, the electrode contact 
window was patterned along the rectangular orbit. Its 
outline is shown by white curves in Fig. [2j For current 
injection, an electrode metal is coated on the cavity top, 
which covers all over the top except for some margin 
of the cavity. The precise shape of the electrode metal 
is defined by r < Mr ((f)) with the margin parameter 
M(< 1). Because of an intervening SiC>2 insulation layer, 
the electrode metal only touches the GaAs contact layer 
through the contact window [5]. With this structure, 
we can selectively excite rectangle-orbit modes, while 
suppressing the excitation of WG-type modes. We 
fabricated samples with M = 0.8,0.85, and 0.95. A 
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FIG. 5: (color online) Far- field emission patterns: experi- 
ment (solid curve) vs wave calculation (dotted curve). The 
experimental data are for a diode with the margin parameter 
M — 0.95. The numerical data are for the rectangle-orbit 
mode shown in Fig. El where rapid oscillation due to inter- 
ference is smeared out. 



scanning electron microscope image of a fabricated laser 
diode with M = 0.95 is shown in Fig. [4{a) 

We tested the laser diodes at 25° C using a pulsed cur- 
rent of 500-ns width at 1-kHz repetition. Figure [5] shows 
a measured far-field pattern for the diode with M = 0.95 
at the injection current 500 mA, together with a wave- 
calculated far-field pattern corresponding to the mode 
shown in Fig. O Both data are normalized so that 
integration becomes unity. Before the normalization of 
the experimental data, uniform background contribution 
due to spontaneous emission is subtracted so that the 
emission level around = coincides with that of the 
wave- calculated data. In Fig. [5j one can confirm excel- 
lent agreement between the experimental and calculated 
data. Far-field patterns for samples with M = 0.8 and 
0.85 are found to be similar to that for M = 0.95, while 
substructures of the two main peaks are different. At 
the current strength 500 mA, we checked from an optical 
spectrum that lasing occurs in multimode. We expect 
that substructures of the main peaks observed in the ex- 
perimental data are due to the effect of multimode lasing 
and thus depend on samples. Sample-specific details will 
be reported elsewhere. 

In order to observe which positions of the cavity light 
is emitted from, we performed a near-field measurement 
of the cavity in lasing operation by using a CCD cam- 
era. The result is shown in Fig. HJb) for a sample with 
M = 0.95, where we can observe scattered light at the 
cavity boundary. In good accordance with the numerical 
results, we can clearly observe that light is emitted from 
two points on each side of the cavity. 

Figure [4](d) shows a near-field photo of a lasing diode 
with M = 0.8. The pattern observed in the margin area 
is considered to be formed by scattered light inside the 
cavity, reflecting the internal modal pattern. We show in 



Fig. 2tc) corresponding numerical data showing the wave 
function of a rectangle-orbit mode superimposed by the 
M = 0.8 electrode metal. Comparing these results, we 
find bright regions of the near-field photo corresponding 
to the left-arm and right-arm parts of the rectangle-orbit 
mode. This further convinces us the excitation of the 
rectangle-orbit modes. 

In conclusion, we have demonstrated chaos-assisted di- 
rectional light emission from the deformed disk semicon- 
ductor microcavity lasers. Theoretical analysis of res- 
onant modes strongly localized along the TIR-confined 
rectangular orbit reveals that their directional emission, 
unexpected from the geometry of the rectangular orbit, 
can be fully explained by tunneling to chaotic orbits even- 
tually escaping from the cavity following unstable man- 
ifolds around the critical line for TIR. In experiments, 
selectively exciting the rectangle-orbit modes, we suc- 
ceeded in observing the chaos-assisted directional light 
emission in both the far and near fields, which provides 
decisive experimental evidence of dynamical tunneling in 
a ray-chaotic microcavity. 
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